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1 Introduction 



The problem of magnetic configurations stability has a long history. In 1600 
W.Gilbert published a treatise "On the Magnet, Magnetic Bodies, and the 
Great Magnet of the Earth" where he proposed that magnets can form the 
noncontact stable systems. 

Since then a number of world known scientists, for example, Newton, 
Earnshaw, Heisenberg, Kapitsa, Braunbeck, Tamm, Ginzburg made their 
contributions to the study of this problem. 

The problem of magnetic equilibrium stability can be naturally divided 
into two tasks of static and dynamic equilibrium. 

Unfortunately, in both of these scientific fields a number of prejudices and 
errors appeared which are not entirely solved until now. 

Concerning the static equilibrium, this resulted in unjustified transference 
of the conclusion from Earnshaw theorem about systems instability in electro- 
statics into the field of magnetic phenomena. Partly this error was eliminated 
in the light of magnetic levitation experiments conducted by Braunbeck and 
Kapitsa- Arkad'ev (i.e. in combination both of magnetic and gravity forces). 

Studies [1] as well as dissertation [2] were devoted to solving the problem 
of static equilibrium of bodies, which interact only via magnetic force. 

Particular prejudices have also penetrated the problem of dynamic sta- 
bility equilibrium in magnetic systems. 

At the dawn of the nuclear age, in the years when research of an atomic 
nucleus has been thriving, magnetic interactions were considered as a possible 
mechanism of keeping particles in the nucleus. In 1941- 1947 Tamm and 
Ginzburg have shown that in the case of two interacting magnetic dipoles 
the orbital motion is impossible, due to the particles falling down the center 
both in classical and in quantum mechanics [3]. In physics this fact was called 
"problem 1/r 3 " and together with Earnshaw theorem, extended to the case 
of magnetostatics, resulted in opinion of "global instability of the magnetic 
systems" for a long time. 

However even after these results famous physicists, J.Schwinger [4] for 
example, had their interest in the magnetic model of matter. 

New splash of interest to the problem of dynamic equilibrium in the mag- 
netic systems resulted in creation by Roy Harrigan a levitron in 1983. In- 
ternet provided wide possibilities for popularization of this unusual toy, and 
also other experiments with magnetic bodies [5]. 
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Against this background the importance of theoretical and experimental 
research conducted by V. Kozorez almost ten years before, in 1974 [6,7] 
seemed to attract unfairly small interest among the world physics community. 

V. Kozorez tried to develope the known idea of the solution of the mag- 
netic systems stability problem based on the consideration that magnetic 
particles are extended presented by Heisenberg in the twenties of last cen- 
tury 

He succeeded in building the experimental prototype where a small mag- 
net accomplished quasiorbital motion up to 6 minutes in duration (this pro- 
totype, in contrast to levitron was not patented). 

As we have already noted in [8,9] his theoretical research had rather 
estimating character, because adequate mathematical apparatus [10-13] for 
studying the stability of such systems has not been developed yet. 

In particular, the condition of stability that he has got for the system 
analogous to the system considered in this article is only one of the three 
sufficient conditions for stability. This condition gives the exact expression 
which reflects the well known fact, that on a considerable distance any mag- 
netic system presents a dipole. 

In respect of the experiment per se, from the philosophical point of view 
an experiment as such or computational modeling in principle cannot prove 
the dynamic stability, but can only give certain reasons in support of the 
stability [14]. 

For the first time the strictly analytical proof of an orbital motion stability 
in magnetic systems is given in [8,9]. Analytical conditions for the stability 
of the system formed by two "magnetic dumbbells" have a complicated form 
and values of parameters for the system's stability area were determined 
numerically. 

Therefore it makes sense to consider a simpler magnetic system that we 
call Orbitron for convenience. Here we analytically prove not only the exis- 
tence of stable orbits but also stability conditions which have simple physical 
meaning for this system. 

In this system a movable body is a small permanent magnet with an axial 
symmetry. Its interaction with magnetic field is described by magnetic dipole 
approximation, and its motion obeys the laws of rigid body motion. 

This differs from the model accepted in work [15]. We do not use the anal- 
ogy originated from the attempts of classical description of such quantum- 
mechanical parameter as particle spin. 

Thus one of tasks of this article is to provide the motion equations of such 
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magnetic rigid body for the systems like Kozorez's prototype and levitron. 

Here, for the mathematical model we call Orbitron, we give accurate an- 
alytical prove of the existence of stable orbital motion of a small magnetized 
rigid body, described as a magnetic dipole based on the theorem from [12]. 
The sufficient conditions of stability in this model have a simple form allowing 
clear physical interpretation. 
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2 Hamiltonian formalism for magnetic dipole 
in the axisymmetrical magnetic field 



Let's consider the Hamiltonian dynamics of "small" rigid body in the ax- 
isymmetrical magnetic field assuming magnetic dipole approximation. Such 
field can be created by cylindrical magnets, solenoids, current-carrying rings 
and other objects with axial symmetry along z axis. 

The variant of such formalism can be obtained from the formalism of 
work [16] by the limiting process m\ — > oo. Thus we consider the first 
body immobile and z axis oriented. Then its dynamical variables disappear 
from consideration or become model parameters. Therefore the group of 
symmetry of the task converges to 5*0 (1). 

For constructing Hamiltonian dynamics based on Poisson structures it is 
necessary to specify a Poisson manifold and also kinetic and potential energy 
of the system. 

Poisson manifold of Orbitron is the direct product of Euclidean spaces 



with Poisson brackets for correspondent generatrix. 

Generatrix for Orbitron will be: Xi - dipole coordinates; pi - its com- 
ponents of momentum (orbital motion); rti - components of dipole intrinsic 
moment of momentum; - components of directing unit vector of dipole's 
axis of symmetry. 

Nonzero Poisson brackets between generatrix on P look like 



Casimir functions of this Poisson structure that are easily checked will be 
V 2 = 1 and (u, n) = const. 

System Hamiltonian we write down in the form: 



(1) 



P = Rl x Rl x Rl x Rl 



(2) 



{rii, Uj} = Eij k v k \ {n h rij} = e ijk n k . 



(3) 



h = T + U(r,c,c,c ), 



where 



(4) 



U = -(fl-B) 
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(5) 



B(f) = B r (r,c')e r + B z (r,c')e z 



(6) 



r = \f\; 
e r = f/\r\; 



e z ■ e r = x 3 /r 
- V ■ e r ; 
= e z ■ V = z/ 3 . 



As usual, kinetic energy of movable body (dipole) consists of kinetic en- 
ergy of both translational and rotational motions [16,8]. 



T(p 2 ,n' 



1 2 

p H — n , 

2AT 2 ' 



where M - dipole mass; a = — ( as we h as before we suppose, that l\ = I 2 = 
J_l, where I\, 13 - intrinsic moments of the body's inertia). 

Get the system of motion equations for magnetic dipole in axisymmetrical 
magnetic field: 



(7) 



< 



' r — p/M; 

f>= -d r Ue r - l -{d,UP1{e z ) + d c „UPi{v))- 
a(n x z7); 
—V x {e r d'i 



v 
n 



e z d c >»)U, 



where P]_ - projector on the plane perpendicular to the vector e r , i.e. P±(e z ) = 
e z — c e r and P]_(v) = v — c e r . 

Expressions for the force and the force momentum acting on a dipole in 
an external magnetic field are well known. We can show that the second and 
fourth equations of the system (7) can be presented in classical representa- 
tion. 

Concerning the second equation in the system (7), it has been obtained 
from the standard expression of Hamiltonian formalism 

(8) ^ = {pi, H} = { Pi , U} = d r U{ Pi , r} + d e ,U{pi, c')} + d c „U{ Pi , c")} 

For potential energy in form (4) we obtain a classic expression of the force 

(8a) $={p,H} = {p,U} = -VU = V(fl-B) 
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Regarding the fourth equation in the system (7), the potential energy 
of a dipole is described by formula (4) in the axisymmetrical magnetic field 
which is described by formula (5), therefore we obtain 

(9) {e r d c n +e z d c ».)U = -fiB 

Then we get the last equation in the system (7) in usual classical repre- 
sentation 

(10) n = i2isxB = jlxB 

Now the system of equations (7) can be written in the form: 

''r = p/M; 

(7a) 

k n = jl x B, 

A few remarks are necessary regarding the systems of motion equations 
(7,7a). 

1. Both systems are correct in the quasi-stationary electromagnetic field 
approximation [17,18]. This approximation is characterized by the possibility 
to neglect the finiteness of electromagnetic disturbances propagation speed 
and displacement current in the range of the system and calculate magnetic 
fields using formulas of magnetostatics. 

2. The system of equations (7) uses the concept of magnetic potential 
energy, which is incident to long-range action conception in classic mechanics. 
As just was mentioned, this is possible in quasi-stationary approximation. 
The chosen form of the potential energy, as in formula (3), describes not 
only dipoles but also wide enough class of axisymmetrical magnetic bodies. 

3. The system (7a) corresponds to the concept of short-range interactions 
in the electromagnetic field theory. Therefore these equations are obviously 
valid not only for the axisymmetrical magnetic field but also describe the 
motion of a dipole in an arbitrary external magnetic field. 

4. We consider a magnetic dipole, as a small magnetized rigid body with 
axial symmetry as in levitron for example. Equation (1) in work [15] in this 
case could not replace the third and fourth equations of the system (7a). 
This distinguishes our mathematical model from that accepted in work [15]. 
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3 Mathematical model of Orbitron 



Not all axisymmetrical magnetic fields can create the possibility for a stable 
orbital motion of a magnetic dipole. For example, the field of magnetic- 
dipole type results in "problem 1/r 3 " mentioned in introduction. Therefore, 
it may be useful to use Heisenberg's hypothesis about the possibility of stable 
magnetic configurations with magnetic extended bodies (see also [7]). 
Here we offer the following model of Orbitron. 

Put two magnetic unlike poles on axis z at points ±h. These poles create 
the axisymmetrical magnetic field in which a magnetic dipole is moving. We 
assume that stable orbital motion of the system is possible under certain 
parameters. 

It is important to give some explanation here. Equations of magnetostat- 
ics do not suppose the existence of isolated magnetic charges. However, the 
field outside a thin solenoid, for example (the same for the thin cylindrical 
magnet) will coincide with high accuracy with the field of two poles [19]. 
On the other hand, the field inside the solenoid not only does not coincide 
with charges field but also opposite in sign, so that the flow through the 
unbounded surface embracing only one pole is equal to zero, as required by 
magnetostatics equations. 

It is assumed that the dipole moves a sufficient distance from the poles of 
the magnet, which is the source of the field, and the model of two magnetic 
charges describes the field with high accuracy. 

Thus, magnetic field in the system has the form of sum of the coulomb 
fields of two charges ±k: 



by formula (5). 

Then for the potential energy of a dipole in the magnetic field we get the 
expression 



(10) 




e=±l 




where each of the fields B £ , and consequently the total field can be presented 



(11) 



U(r,c',c",c"') 




eU E (r,c ',c ",c '"), A = noKfi 

E=±l 
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where 

// . /// 

12 U E (r,c ,c ,c =^7 r^- 

R £ (r,c ) d 

and 

, / , \ 1/2 

(13) R £ (r, c ) = (r 2 - lehrc + h 2 ) 

Let's show the first derivatives of the function U £ : 

c" 3 (re — ehc") (r — ehc) 



(14) d r U £ 



(15) d,U t 



R £ (r,c') 3 R £ (r,c'f 
3 (re — ehc") ehr 



Re(r,c'Y 



(16) djU. 



R £ (r,c ') 3 



(17) d,„U £ = 



RAr, c') 3 
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4 Example of stable orbital motion 



The main aim of this work (i.e. Part I) is to prove the existence of stable 
orbital motion in the systems of bodies, which interact only by magnetic 
forces. The example of such system is described in a section 3, and example 
of a stable orbit will be the circular orbit in plane z — 0. 



4.1 Relative equilibrium 

A special role in orbital motions stability of Hamiltonian systems plays the 
so-called relative equilibrium [11,12], i.e. such trajectories of the dynamic 
system which simultaneously are one-parameter sub-groups of the system's 
invariance group. 

As it has been already mentioned, the invariance group of Orbitron is 
SO(l). Every one-parameter sub-group of this group is characterized by the 
intrinsic rotational angular velocity u = ue z . For our problem the rate of 
change of any physical value v along the orbit of the sub-group will be set 
by the formula v = Co x v. 

Therefore, for the relative equilibrium to exist the following relationships 
must hold 

' r — uj(e z x r); 
p = u(e z x p); 
V = u(e z x z7); 
= u(e z x n). 

We show that a dynamic orbit for which these relationships are satisfied 
exists. Examine an orbit, spatially located in the z = plane. Also suppose 
that V || e z and n || e z . Then c=c = 0,c =±1 along the whole trajectory 
and d c 'U = d c "U = d c »>U = as follows from formulas (11,15-17). 

So, the third and the fourth equations of the system (7) then hold iden- 
tically, the first and the second are reduced to the second order equation: 

(19) MF+f^) r = 

V ' / \r=ro 

On condition that {d r U)\ r=ro > equation (19) has solution corresponding 
to the motion on the circumference with radius r and frequency, which is 



< 
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determined by relationship 



(20) 



( 



d r U 



r 



) 



u 2 M 



Thus, one can prove that the reduced orbit indeed is a relative equilib- 
rium. 

Theorem 4.8. in [12] is a suitable instrument for investigating stability of 
relative equilibria on Poisson manifolds. Important advantage of group the- 
oretical methods is that the functional space of investigation of trajectories 
is substituted by investigation of finite-dimensional vector space of dynamic 
variables variations in a fixed point on the trajectory. Thus the investigation 
approach for stability is very similar to the study of a function's conditional 
extremum by Lagrange multiplier method. 

4.2 Choice of supporting point 

Lets set the point on an orbit of relative equilibrium 



Notice that we do not fix the sign of the mechanical moment n , it can 
be arbitrary. As for a sign of po, for a positive angular velocity its value will 
be positive. 

Lets show that in supporting point the following relationships hold 



(21) 



'x = r ei; 
Po = Poe 2 ; 
v = -e 3 ; 



{n = n e 3 ; 



(22) 




0; 



d c "U\z. = 0; 



Since in supporting point 



(23) 




therefore 



(24) 
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From the expressions of potential energy derivatives (15-16) 

(25) (d c 'U £ ) [ze = + 2 , {d c nU £ )\ Ze = ^ 2 - 2 

notice that both expressions do not depend on e, meaning that in sum on e 
(with e- multiplier) they will give 0. 
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4.3 Necessary condition of stability and Lagrangian 
coefficients 

As motion integrals we will take 

{ja = %iP2 ~ X2P1 + n 3 ; 
C 2 = \ 2 {p,n); 

where 1st line represents a third conserved quantity of a body total angular 
momentum, and the other two are Casimir functions of the system. 
Write out the correspondent differentials in z e point 

{{dh)\ Ze = Podx! + r dp 2 + dn 3 ; 
{dC x )\ Ze = -A lC ^ 3 ; 
(dC 2 )\ Ze = \ 2 {n dv 3 - dn 3 ); 

Efficiency function (adjoined Hamiltonian) looks like 

(28) H = T + U -uj 3 + Xxd + \ 2 C 2 

The necessary condition of stability in theorem 4.8. [12] requires the 
differential of efficiency function to be equal to zero in a supporting point, 
i.e. dH\ = 0. 

|ze 

For the differential of potential energy we have 

(29) dU\ Ze = d r Udxx + d d n\Jdv 3 
For the differential of kinetic energy we have 

(30) dT\ Ze = j^dp 2 + an dn 3 ; 
Collecting the differentials of efficiency function, we get 

(31) dH\ Ze = {d r U\ Ze - ujp )dx 1 + - ur^j dp 2 
+(d c "'U\ Ze — Ai + \ 2 n )dv 3 + (an — uj — \ 2 )dn 3 
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Equating dH\ Ze = 0, we derive the following expression for Lagrange 
multipliers 



(32) 

where 
(33) 



'po/M = wr ; 
cupo = d r U w = 
X2 = ano — lo; 

Ai = d c mU\ Zs + A 2 n = K + n (an - u), 



K = dr»>U\. 



2irR 3 



the first equation in (32) is an ordinary relationship between linear and 
angular velocity during circular orbital motion. 

second equation in (32) represents the equality of centrifugal (on the left) 
and centripetal (on the right) forces. 

From this two expressions we get the relationship for angular velocity, 
namely: 



(34) 



1 IK 

Mu 2 = -d r U lZe = ^ 
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4.4 Allowable variations 



For the application of the sufficient condition of stability in the theorem 4.8. 
in [12] it is necessary to extract a linear subspace of allowable variations. 

Let's consider the variations of the dynamic variables annihilating the 
differentials in formula (27). 

From the second line in (27) it follows, that 5u 3 = 0, then it ensues from 
the third line, that 5n 3 = 0. 



Thus, we obtain 



(35) 




0; 
0; 



Hence it ensues that the variations in the form 
(36) Sx 1 , 5x 2 , 5x 3 ; Spi,Sp 3 ; 5v 1 ,8v 2 \ 



Sni, Sn 2 



can be considered as independent variations, furthermore, we must exclude 
from this subspace the direction which is tangent to the orbit 

It ensues from formula (18), that this direction (in z e point) is determined 

as 



(37) 



In order to eliminate the variation (37), we impose another additional 
condition on variations, and then we get the constraints 

(5u 3 = 0; 

Sn 3 = 0; 

Spi = g&fc; 
l$P2 = 



' 5x = 


r e 2 ; 


Sp = 




* 8V = 


0; 


^Sn = 


0. 



(38) 



and an independent set of variations will be 
(39) 5x\5x 2 ,5x 3 ; 5p 3 ; 8v\5v 2 ; 



8ni, <5n 2 ; 
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4.5 Basic quadratic form 



Sufficient condition for a minimum consists in positive defmiteness of quadratic 
form of type d Hi (6z, Sz ), where variation vectors Sz, 5z must be expressed 
through independent variations (39) taking into account the constraints (38). 
Quadratic form defined in independent variations we denote by Q. 

Calculations of the efficiency function hessian (adjoined Hamiltonian) and 
basic quadratic form in independent variations were performed in Maple. 

For better structuring of the expressions indefinite Lagrange multipliers 
are hidden at the first stage. 

After insignificant transposition of columns (and corresponding lines with 
the same number) the matrix of basic quadratic form acquires a form 



(40) 



Qn 


























Q22 


























Q44 


























Q33 


Q35 




















Q35 


Q55 


Q57 




















Q57 


Qn 




























Ch 






















Q s 



Lets write out non zero elements from matrix of quadratic form (per line) 



(41) 



Q 
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h 2 - Ar 2 Q K 
W 2 B 2 



+ Mou- 



rn 



Q22 = 3 (J^ + Muj2 



(42) 



Qa4 ~ M 



(43) 
(44) 



Q 



38 



3r 2 - 2h 2 3K 
B 2 'B 2 ' 



Q 



35 



^55 — ^1; ^57 — ^2 



R 2 
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(45) 



Q 77 = a 



(46) Q m = Ai = Q 55 , Qes = A 2 = Q 



57 



(47) Q 88 = a = Q 77 

Substituting Mu 2 for expression (34) in Qn,Q22, we obtain 

,a, \ „ 3K4h 2 -r 2 

(41a) Q« = 

12K 



(42a) Q 



22 



R 2 
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4.6 Conditions of positive definiteness of the basic quadratic 
form 



For matrix Q to be positive definite it is foremost necessary that all diagonal 
elements of the matrix are positive. Q22, Qaa, Q77, Qss are scienter positive. 
Remaining conditions are as follows 



(48) 




~r? W 

3K 3r 2 ~2h 2 _ 
R 2 R 2 ' 



= Qm = Ai = K + n (cm - w); 
The first two conditions result in purely geometrical limitations 



(49) 



(Qn > 0) & (Q 33 > 0) 




3 h ] \h 



These conditions of positive definiteness of the matrix Q have to be sup- 
plemented now by the conditions of positive definiteness of two submatrices 
of 3 x 3 and 2x2, namely 



(50) 

and 
(51) 



Q33 Q35 
Q35 Q55 Q57 
Q57 Q77 



Qm Q68 

Q&8 Q8S 



Taking into account the above mentioned considerations it is sufficient 
for the matrix (51) to check the condition of positiveness of its determinant 
of QmQ 88 ~ Qls- Thus, additionally to the conditions (49) the following 
condition is added 

(52) < Q m Q 8 8 ~ Qls = aK + u(an - u) 

Now we investigate the conditions of positive definiteness of the matrix 
(50). The first condition Q33 > we have considered already. 
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Thus the additional conditions of positive definiteness of matrix (50) are 
reduced to positiveness of two determinants 

(53) < Q 33 Q5,5 - Q35 
and 

(54) < Q33Q55Q77 - Q33Q57 " 
Condition (54) can be also written in form 

(54a) Qrr(QaaQsB - Q35) > Q33QI7 

and, since Q77 > then (54) transforms into (54b) which replaces condition 
(53) as it accounts for it 



(546) 



Q33Q 



55 



Q 



77 



57 



So, condition (53) is superfluous and it is necessary to study only condi- 
tion (54). Remind that 



(55) 



33 



3K 3rg-2/t 2 . 
R 2 R 2 ' 



Q 

n — o Krp . 

V35 — R 2 ' 

Q55 = K + n (an - u) = Q 66 ; 
Q 57 = an -uj = Q 68 ] 
Q 77 = a = Qss] 



Now write down (54) in form 
(54c) Q55Q77 - Q57 > 



Q77Q35 



Q 



33 



as supposed Q 33 > 0. 

Taking into account formulas (55), condition (54c) is equivalent to 



(54d) 



QmQ88 ~ Qls > 



Q77QI5 

Q 



> 



33 



Thus, condition (52) is a consequence of condition (54) and Q33 > from 
(48). It means that condition (52) can be omitted. 
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We have the following reduced number of conditions for matrix Q positive 
definiteness: 



(56) 



U <- (^11 - fl2 fi 2 , 

< Q55 = Qee = Ai = K + n (cm - w); 
< Q33Q55Q77 — Q33Q57 ~ Q77Q35 



We investigate condition (54) in form 

< Q 3 3(Q5 5 Q77 - Q57) 

= Q33(Q66Q88 — Qls) 



Q77Q35 



Q77Q 



35 



^[(3rg - 2/j>(an - u) - 2aKh 2 } 



That is 



a 



an -u>)> K^ 2 



2h 2 



2h 2 



So, condition (54) is equivalent 



(57) 



K 



-(ano-u) > 3 



2\ h> 1 



In particular, as Q33 > 0, we have ario — 00 > 0, and it means that 
conditions Q55 = Qm — Ai > are fulfilled a priory and can be omitted. 
Therefore, conditions 



(5* 



2 < a < 9- 

3 ^ ft ^ z ' 



K 



define positive definiteness of form Q. 
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4.7 Physical meaning of positive definiteness condi- 



In the previous section we have established the conditions of the system's pa- 
rameters which provide positive definiteness of basic quadratic form, namely: 



The first condition in (59) is purely geometrical and determines a possible 
range for a radius of the orbit representing a relative equilibrium (21). 

The second condition is dynamic and determines lower boundary for the 
intrinsic moment of momentum of the body. In particular, it means that a 
body must be sufficiently rapidly revolved. Therefore it is worth to solve this 
inequality relative to n . 

Using relationships (34), we obtain 



Value of - corresponds to the intrinsic moment which a body would have 
if it were to revolve with angular velocity u athwart to the own axis of 
symmetry. 

Value of u(MrQ) is simply the orbital moment of momentum (L z )\ Ze . 

For condition (60) one can give such physical meaning: intrinsic moment 
of body rotation must be of the same or higher order of its orbital moment. 

Indeed, multiplier before u(MrQ) is a geometrical factor which is ^ = 1.0 
equal to ~ 1.33, and at y = 1.5 equal to ~ 0.2. 

In these estimations the first term in the right part of expression (60) can 
be neglected. 



tions 



(59) 




(60) 
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5 Numeral simulation 



From a mathematical point of view the stability conditions obtained here 
allow a wide range of parameters values of the problem to exist, however not 
all of them can be physically realized. From a physical perspective the values 
of parameters are limited by the properties of present materials. 

In addition it seems difficult to realize in practice high rate of rotations, 
especially as far as it concerns intrinsic angular velocity of movable magnetic 
body (dipole). 

Therefore it appears necessary to specify such values of parameters which 
can be realized in an experiment. 

For the magnets, made from Nd—Fe—B, we have the following character- 
istics: p = 7A-10 3 (kg/m 3 ) - density and B r = 0.25(T) - remaining; induction. 
Then it is easy to obtain magnetic "charge" of the poles k = 17. 6 (A • m). 
Distance between the poles is L = 2h = 0.1(m). 

A movable magnet we choose in a form of cylinder (disk) with the diam- 
eter of d = 0.014(m) and height / = 0.006(m). Then disk magnetic moment 
H = 0.1801 -m 2 ). 

As a result for the orbit with the radius = 1.5h = 0.075(m) we obtain 
the angular velocity of the orbital motion u = 1.54(rad/ sec), with minimum 
angular velocity of disk intrinsic rotation in this case is Q = 72.8(rad/sec). 
Such values of angular velocity appear fully reasonable. 

Using the indicated values of Orbitron parameters the numeral modeling 
of orbital motion was conducted under the deviations of initial values of 
dynamic variables from the values, which correspond to relative equilibrium 
within 1% error. 1000 castings which showed the stability of orbital motion 
were accomplished by Monte Carlo method (i.e. by random selection of the 
initial values in the vinicity). We will elucidate this in more detail in the 
next parts of this work. 
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6 Summary 



The main aim of this work was to give constructive proof of stable orbital 
motions existence in the systems of bodies, which interact only by magnetic 
forces. 

For this purpose it is enough to analytically prove the existence of stability 
for one orbit in comparatively simple system described by equations which 
do not contradict the laws of electrodynamics and classical mechanics. 

We named such a system Orbitron, found its parameters which can be 
physically realized and conducted the Monte Carlo numeral modeling. 

To be continued . . . 
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